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It is expected that, after a quench, an interacting quantum system eventually reaches thermal 
equilibrium. Here we show that this is not always the case. We investigate the dynamics of the 
order parameter of a strongly coupled superconductor after a quench by holographic techniques. The 
gravity dual that we employ is the AdSs Soliton background at zero temperature. Time evolution is 
first investigated by computing the quasi normal modes (QNM) associated to the superconducting 
order parameter. QNM's are purely real which suggests undamped time oscillations of the order 
parameter which prevents thermalization of the dual field theory. Results for the time evolution of 
the order parameter after an abrupt change of the chemical potential are consistent with this picture. 
By tuning the quench strength we identify a region of parameters for which the oscillating pattern 
of the order parameter becomes increasingly intricate. We also discuss the range of applicability of 
these results and the physical origin of the lack of thermalization. 
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Advances in cold atom physics and numerical tech- 
niques have renewed the interest in the route and con- 
ditions for thermalization in closed, strongly interacting, 
quantum systems after a quench. Until recently the- 
oretical progress had been relatively slow as traditional 
analytical and numerical approaches are not in general 
well suited to describe far from equilibrium dynamics. 
However the application of the holographic principle, also 
called the AdS/CFT correspondence has opened [3- 

new research avenues to tackle this problem. 
This correspondence states that, under certain condi- 
tions, strongly coupled gauge conformal field theories 
(CFT) in d dimensions are dual to a classical gravity 
theory in Anti de Sitter (AdS) space in (i+ 1 dimensions. 
The out of equilibrium dynamics of the field theory has 
an especially appealing interpretation in the gravity dual: 
it corresponds with the time evolution of a mass distribu- 
tion in an asymptotic AdS background. Thermalization 
is thus related to the formation of a black hole in the 
gravity dual 0, B- In the context of holographic su- 
perconductivity [I4I there are recent studies on the time 
evolution of the order parameter 0, [^, [H, [T^l at finite 
temperature. In [T^ it was found that, in the probe limit, 
QNM of the order parameter have finite real and imagi- 
nary parts which indicates that, for long times, the super- 
conducting order parameter oscillates with an amplitude 
that decays exponentially as thermal equilibrium is ap- 
proached. Moreover Ref. [lj| also identified the Goldstone 
mode that signals the spontaneous breaking of the U{1) 
symmetry at sufficiently low temperature. These results 
have been recently confirmed in [!|] by the explicit calcu- 
lation of the order parameter after a quench taking into 



account backreaction effects. Similarly a previous study 
Q, also including backreaction, showed that, for suffi- 
ciently long times after a quench, deviations of the order 
parameter from the thermal equilibrium prediction are 
exponentially small. Moreover a quench in a Reissner- 
Nordstrom-AdS black hole can induce an instability in 
the metric that leads to holographic superconductivity 

i- 

Recent experimental and theoretical [l|, [l5|, [l^ re- 
sults in condensed matter points to a richer phenomenol- 
ogy. Integrability and localization P, [3, [l^ can prevent, 
or slow down [15| . thermalization after a quench at zero 
temperature. In the context of weakly coupled supercon- 
ductors [1^ undamped oscillations of the order parame- 
ter after a quench have been observed provided that the 
final coupling is much larger than the initial one. For 
intermediate quenches the superconducting gap is oscil- 
latory in time with an amplitude that decays towards 
equilibrium in a power-law fashion [l6j . 

Motivated by these results we explore novel dynami- 
cal regions in the time evolution of holographic super- 
conductors. Our main aim is to identify a range of pa- 
rameters for which thermalization does not occur. For 
this purpose we employ a superconducting AdS Soliton 
background [13, [ill which has a well defined limit at zero 
temperature. This is a key departure from previous holo- 
graphic studies [1, [^, [l3| which were carried out at finite 
temperature, namely, in an AdS black hole background. 
At finite temperature thermalization always occurs af- 
ter sufficiently long times. Therefore it is important to 
work in a background which has a well defined limit at 
zero temperature. In order to gain insight into the late 
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FIG. 1. (Color online) First (blue) and second (green) order 
QNM Refoj] © as a function of ^i. A vanishing QNM (black 
dot) signals the superconducting transition at /j. — fic ^ 1.88 
[l7| . Within the numerical precision of the calculation the 
imaginary part of the QNM vanishes for any fi. This suggests 
that after a quench the dual field theory does not thermalize. 
The dynamics of the order parameter (see Fig. [3]) becomes 
more intricate for fi > 3.8 which corresponds to the minimum 
fi needed to excite the second order QNM. 



time dynamics of the order parameter we study first the 
QNM [1^ related to the order parameter [2^. We have 
fomid that in the AdS Solition background the QNM's 
are always real which suggests that the order parameter 
oscillates without any damping. A detailed analysis of 
the time evolution after a quench confirms this predic- 
tion. Therefore no thermalization occurs, namely, after 
a quench the strongly coupled superconductor does not 
reach thermal equilibrium. However the oscillating pat- 
tern strongly depends on the quench details. In some 
cases we observe quasi chaotic oscillations. 
The model. - 

We study the time evolution of a holographic supercon- 
ductor in an AdS5 Soliton gravity background. This 
background [3| is constructed from the usual AdS5 ge- 
ometry by a double Wick rotation followed by compacti- 
fication of one of the spatial dimensions. In order to have 
a smooth geometry it is necessary to impose periodicity 
X X + 7~ in. the compactified dimension. The resulting 
geometry has no horizon, it looks like a cigar with a tip 
at r = tq. The most stable configuration which satisfies 
this constraint, the so called AdS Soliton [31 , is fully 
specified by the metric. 



dx^+dy^)+f{r)dx^ (1) 



with f{r) = — p-. In this background we in- 
troduce a charged scalar field ^ minimally coupled 
to a Maxwell field, which leads to the following five- 
dimensional Einstein-Maxwell-scalar gravity theory. 
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We refer to ([2|) , together with the metric ([T]) , as the AdS 
soliton holographic superconductor We note that 

this model, unlike the usual AdS-Schwarzschild super- 
conductor, has a well defined zero temperature limit. 

The properties of this model, investigated in in the 
probe limit and in |2l| including the backreaction of the 
scalar on the metric, depend on the value of the chemical 
potential fi, which corresponds to the time component of 
the gauge field at the boundary (see below). For fi < fic 
the conductivity in the linear response regime vanishes 
and therefore the field theory resembles that of a Mott 
insulator (T^]. However for fj. > fic the scalar condenses 
and the system becomes a superconductor though the 
background is still an AdS Soliton. It is therefore ex- 
pected that some of the Mott insulator physics might still 
be at play in this regime. For sufficiently large ^ /i^ a 
confinement-dcconfincmcnt transition occurs for any fi- 
nite temperature but not at strictly zero temperature 
[2l| . Superconductivity survives this transition but the 
gravity background becomes an AdS black hole [l2| . 
We restrict ourselves to the probe limit in which the back- 
reaction of the gauge field and scalar on the metric is 
negligible. Several reasons justify this choice. It was re- 
cently shown in [2^ that the insulator-superconductor 
transition in the AdS Soliton does not proceed from a 
dynamically instability. Therefore it is not in principle 
possible, even in a calculation including backreaction, to 
induce the transition by a quench. Moreover the results 
of including backreaction, at zero temperature are 
similar to those obtained in the probe limit [l7| . 
The equations of motion (EOM).- 

Since we are interested in time-evolution we seek solu- 
tions of the EOM that depend on both the temporal 
and radial coordinates: A = (^t, A^, 0, 0, 0), * = |^|e'^ 
where and 9 are functions of t and r [23j . 

In the following, we will work with the gauge-invariant 
quantities AI = A — dO. Without losing generality we 
set L = tq = q = 1 and — —15/4 so that the time 
evolution depends only on the chemical potential. For 
numerical purposes it is more convenient to make the 
following change of variables, z = 1/r and \^\ = |^'|/r'^/^ 
so that the boundary is at z = 0. The resulting EOM 
are: 



+ {z^-l) z^dM + ^'dfli^l 

- 1) zdtd^Ad 
+ {z*-l) z^dlMt + 2z^Mt\ij\^, 



(3) 

= (z^ - 1) zdtd^M, + (z^ -f 3) dtM, + (3z^ + z^) d,Mt 



^d^M, + 2zA4|^|2 + z^dtd.Mt, 



(4) 
(5) 



m^l^-P) 



= - (z** - 1) z\iP\d^M^ + z^\i:\dtMt + 2z^Mtdt\ip\ 

- 2M, (^z (z^ - 1) 9,1^1 -f (z^ + 1) . (6) 

We note that Eqs. ([1]), (O and © arc not independent 
since Eq.® can be obtained from Eqs. @ and 
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Boundary conditions. - 

In order to solve the EOM above it is necessary to specify 
initial conditions (t = 0) and boundary conditions at the 
tip z = 1 and at the boundary 2 = 0. Near the boundary 
and with the definitions above, 

1^1 = \Mt)\ + \Mt)\^ + o{z^), 

= (7) 

where ^ and p stand for the chemical potential and the 
charge density of the dual field theory respectively. We 
choose solutions for which |'0i| = 0, Itjj^l ^ 0. Following 
the standard AdS/CFT dictionary [12| the condensation 
of the operator O2, 



mt)) = \Mt)\ 



(8) 



stands for the order parameter of the dual field the- 
ory which describes the time evolution of the strongly- 
coupled superconductor. Both the gauge field and the 
scalar must not be singular at the tip. A simple choice is 
\4>\ — ai + 02(1 — z) + . . . , Mf =03-1- 04(1 — z) + . . . and 
Mz = 05 + 05(1 — z) + . . . , where the a^'s depend only 
on time. In order to obtain the initial conditions, i = 0, 
we solve the EOM in the static limit for p,{t ~ 0) = pi. 
QNM of the scalar field. - 
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FIG. 2. (Color online) (C'2(t)> (O for different quenches, 
defined hy p. f,jj.i > fic- In full agreement with the QNM 
results (see Fig. [1} only undamped oscillations are ob- 
served, (a) Pi — S,pf = 3.2 (red), pi = 2.6, pf = 2.7 
(green), pt = 1.89, /i/ = 2.0 (black). The frequency of 
the oscillations is controlled by pi. (b) pi — 1.89 and 
Pf = 1.90(blue),2.0(black), and 2.3(dotted). The ampli- 
tude of {02{t)) increases with \pf — pi\. (c) Blown up of 
Pi = 1.89,^/ = 1.90 in (b). The oscillations are small but 
clearly visible. These results indicate that there exist strongly 
interacting field theories that after a quench do not thermal- 
ize. 

In order to gain insight on the asymptotic dynamics of 
(02(0) compute the QNM associated to the scalar 
1-01. The QNM's in the gravity theory correspond to the 
poles of the retarded Green's function of the dual field 
theory [l^. Therefore it is a powerful tool to investi- 
gate the asymptotic time evolution of the order param- 
eter after a small perturbation. We define = 



{tpoiz) + e~*'^*'i/'(z)| where ijj is small with respect to the 
static solution tpo{z). We carry out a similar expansion 
for the other fields. We keep only terms in the EOM that 
are leading in these small perturbations. The QNM are 
defined as the discrete set of complex. 



uj{p) 



(9) 



frequencies that solve the EOM. A zero QNM is a 
signature of a phase transition [25j . As was expected, 
see Fig. [TJ it occurs at p = pc where the insulator- 
superconductor transition occurs [TtI . 
For p < pc (region I) , the EOMs decouple and the 
relevant QNMs correspond to the pole of the retarded 
Green's function of the scalar. As shown in Fig. [1] the 
QNM's only contain a real part with ujr oc \p — pc\. 
For p > Pc in the superconducting phase, the calcula- 
tion is more challenging since the scalar is coupled to 
a gauge field. In order to compute the QNM's we use 
Chebyshev spectral method where the perturbations 
above are explicitly gauge-invariant. This is important 
in order to avoid spurious QNM's. In Fig. [T]we plot the 
two lowest order QNM's w as a function of p,. The black 
dot is the marginal stable mode which corresponds to 
cj = 0; this is the critical point. For all p > pc, cui = 




FIG. 3. (Color online) (02 (t)) dS]) for pt = 5.0 and, from 
top to bottom, Pf = 5.04,5.07,5.14. For any pi > 3.8 and 
sufficiently small quenches oscillations of {02{t)) become in- 
creasingly chaotic. For smaller p's the second order QNM 
cannot be excited and oscillations are not chaotic. 

and UJR 7^ which indicates undamped oscillations of 
{02{t)). This is a striking indication that a holographic 
superconductor in a AdS Soliton background does not 
thermalize. As can be observed in Fig. [T] there is a 
jump of the QNM at /i — pc- We believe that this 
reflects the effect in the superconducting phase of the 
mass gap. which according to [l^, is related to Mott 
physics. Region II is also characterized by a growing 
(decreasing) lowest (second) QNM's as p increases. This 
trend is reversed for ^ > 3.8 (region III). In the next 
section we shall see that precisely around these values 
there is a qualitative change in the oscillating pattern 
of the order parameter, from simple sinusoidal to more 
complex oscillations. In order to fully confirm this lack 
of thermalization we study {02{t)) next. 
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Full time evolution of {O2) after a quench. - 
QNM's only provide information about the dynamics for 
asymptotically long times and for small perturbations. 
In order to study shorter times and stronger pertur- 
bations we compute explicitly the full time evolution 
of {02{t)) (IH) after a quench of the chemical potential 
ji{t) at t w 0. This type of quench has a clear physical 
interpretation: on the field theory side, a larger fi 
indicates a stronger order parameter [l^. A larger 
amplitude of {02{t)) indicates a stronger interaction. 
Therefore by changing abruptly the chemical potential 
we are also modifying the effective interaction. At least 
qualitatively this is exactly the type of quench which is 
commonly used in condensed matter [l6j . 
For numerical reasons the quench cannot be instan- 
taneous. We assume that goes from /i(0) = 
to fi{T) w fif in a time r — 0.5/ fii (time is measured 
in units of fii) which is much smaller than the relax- 
ation time of the system. The exact form of the time 
dependence for < t < r is not important. For the 
sake of simplicity the quench is defined by imposing an 
additional boundary condition at //(t) ~ fJ-f. The latter 
is not strictly a quench as the effective value of the 
interaction strength changes smoothly on time in order 
to verify the boundary condition at t = t however it is a 
good approximation since t ~ 0.5 tf ^ 20. 
The EOM (g]), © and © with boundary conditions © 
and the quench definition above were solved by using 
spectral methods. We discrctize the EOM on a two 
dimensional Chcbyshev grid with 100 points along the 
t direction and 20 points along z direction. We study 
times up to = 20. Technically it is feasible to go up to 
tf = 40 by adding more points along the t direction but 
the results are very much the same as those for = 20. 
In Fig. Hand Fig. [3] we plot (02(0) ® for different 
quenches strength parametrized by fj,{t = 0) = /i,; and 
fi{t = t) = fif . We explore a broad range of parameters 
with the only restriction that fJ.i,fif > jJLc- Several 
comments are in order: a) undamped simple sinusoidal 
oscillations are observed for jii < 3.8, b) the frequencies 
of the oscillations increase with c) the amplitude is 
an increasing function of j/x,; — This is expected as 
a stronger quench indicates a larger difference between 
the initial and the final coupling, d) similar results are 
obtained by different types of quenches and different 
choices of r, e) similar results are also observed for 
weakly coupled superconductors [l^ in a certain region 
of parameters, f) for fii > 3.8 and small quenches we 
observe (see Fig. [3]) a more complex oscillating pattern. 
As the quench strength increases the simple sinusoidal 



pattern corresponding to fj,i < 3.8 is gradually recovered. 
This minimum fi ^ 3.8 is the one for which the lowest 
two order QNM's contribute simultaneously to the order 
parameter dynamics. In summary thermalization never 
occurs in our strongly coupled superconductor. This is 
the main result of the paper. 

As was mentioned previously the probe limit is enough 
because the transition is not induced by a dynamical 
instability |22| and because previous results including 
backreaction [2l| were qualitatively similar. As a further 
verification we have reproduced in the probe limit the 
main findings of [oj in the broken phase of the AdS black 
hole background which included the backreaction of the 
scalar on the metric. Finally we stress that our results 
are valid provided that quantum gravity corrections are 
neglected. For sufficiently long times it is expected that 
these corrections will induce some sort of thermalization 
in the system. However, even if thermalization finally 
occurs, its typical equilibration time is much larger than 
in other field theories with a gravity dual H, Q . 

In conclusion, we have investigated the time evolution 
of a holographic superconductor by computing the 
QNM's and the explicit time evolution of the order 
parameter after a quench. Our main result is that 
the order parameter does not thcrmalizc. It performs 
undamped sinusoidal oscillations typical of insulators or 
integrable systems which do not reach thermal equilib- 
rium after a quench. We have also identified a region, 
related to the existence of relatively close QNM's, for 
which the oscillations of the order parameter follow 
an intricate pattern. Since the dual field theory has a 
mass gap we speculate that the physical mechanism that 
prevents thermalization is related to localization caused 
by interactions. 

Note: At the time we were finalizing this work a paper 
[26| was posted on arXiv which also studies a quench in 
an AdS Soliton background. 
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